Image-charge induced localization of molecular orbitals at metal-molecule interfaces: 

Self-consistent GW calculations 
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Quasiparticle (QP) wave functions, also known as Dyson orbitals, extend the concept of single- 
particle states to interacting electron systems. Here we employ many-body perturbation theory in 
the GW approximation to calculate the QP wave functions for a semi-empirical model describing 
a 7r-conjugated molecular wire in contact with a metal surface. We find that image charge effects 
pull the frontier molecular orbitals toward the metal surface while orbitals with higher or lower 
energy are pushed away. This affects both the size of the energetic image charge shifts and the 
coupling of the individual orbitals to the metal substrate. Full diagonalization of the QP equation 
and, to some extent, self-consistency in the GW self-energy, is important to describe the effect 
which is not captured by standard density functional theory or Hartree-Fock. These results should 
be important for the understanding and theoretical modeling of electron transport across metal- 
molecule interfaces. 

PACS numbers: 71.10.-w,73.20.-r,72.10-d 



I. INTRODUCTION 

The independent-particle approximation and the asso- 
ciated one-electron orbital picture forms the basis of our 
understanding of chemical bonding and electronic energy 
levels in solids and molecules. The most widely used 
approximations of this type are Hartree-Fock (HF) and 
density functional theory (DFT)ii2. Although the single- 
particle orbitals derived from such schemes do not have 
physical meaning, apart from the fact that the exact DFT 
orbitals generate the exact groundstate density, they are 
routinely used to calculate and interpret physical quan- 
tities of various types. In strongly correlated systems 
such an approach clearly breaks down. However, even in 
weakly correlated systems where the single-particle pic- 
ture is valid, there is no guarantee that the orbitals gen- 
erated by the standard one-electron schemes are those 
which best resembles the true many-body excitations. 

Quasiparticle (QP) wave functions provide a rigorous 
generalization of the concept of single-particle orbitals to 
interacting electron systems. The QP states and energies 
are solutions of the QP equation 3 - 



(1) 



Here Hq is the non-interacting part of the Hamiltonian 
including the Hartree field while T, xc is the non-local and 
energy-dependent exchange-correlation (xc) self-energy 
operator^. The QP energies represent the possible en- 
ergies of a particle (electron or hole) added to the TV- 
particle groundstate, and the QP wave function describes 
the probability amplitude for finding the added particle 
at a given position. (A precise definition and interpreta- 
tion of the QP energies and wave functions will be given 
later in this paper.) 

The GW approximation 3 to Y, xc (both in its self- 
consistent and non self-consistent form) has been 
succesfully used to calculate QP energies of solids^—, 



molecule s 11 ' 12 , and more recently solid-molecule 
interfaces^ 3 -"— . The latter class of systems is par- 
ticularly challenging to describe due to its highly 
inhomogeneous nature: Since the QP states describe 
the charged excitations, the QP energies of an adsorbed 
molecule are strongly affected by the metal surface 
through long range polarization effects (image charge 
effects) which decay as 1 / z with z being the distance to 
the surface, see Fig. 1(a). The inability of any available 
DFT functional to account for this renormalization of 
molecular energy levels, reflects the highly non-local 
nature of the phenomenon. We stress that this does not 
imply that the single-particle picture is invalid in such 
cases; only that the correct (QP) orbitals and energies 
cannot be obtained from a simple mean field potential. 
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FIG. 1: (Color online) Schematic illustration of the change in 
the frontier QP energies and wave functions of a molecule ap- 
proaching a metal surface. The closing of the HOMO-LUMO 
gap due to image charge screening is associated with a change 
in the shape of the orbital. 

Most applications of the QP equation have focused on 
the QP energies while the QP wave functions have been 
much less studied. In fact it is very often assumed that 
the latter, apart from normalization, are identical to the 
orbitals obtained from DFT (V>u)- Under that assump- 
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tion, QP energies can be obtained from first order per- 
tubation theory involving only the diagonal matrix ele- 
ments (V^l^-xc^)!^) thus greatly reducing the compu- 
tational cost of solving the QP equation. 

As an example where the shape of the QP wave func- 
tions play a key role we consider the case of electron 
transport through a molecule connected to metallic elec- 
trodes. The conductance of the junction depends mainly 
on two factors, namely (i) the position of the molecule's 
frontier (QP) energy levels relative to the metal Fermi en- 
ergy and (ii) the overlap between the molecule's frontier 
(QP) orbitals and the extended states in the electrodes. 
The importance of (i) has been studied in detail using 
scissors operator techniques to correct the DFT energy 
levels while keeping the DFT orbitals fixe d 18 i 19 . In con- 
trast, the question of how well the DFT orbitals resem- 
ble the true QP orbitals and the consequences for charge 
transport has only been indirectly studied^— . 

In this paper we show, using many-body perturbation 
theory in the GW approximation, that the QP wave func- 
tions of a molecular wire in contact with a metal surface 
can be qualitatively different from those obtained from 
an independent-particle approximation. While the latter 
method yields orbitals that remain delocalized over the 
molecule upon coupling to the surface, our GW calcu- 
lations show that image charge effects not only causes a 
reduction of the QP energies as previously demonstrated, 
but also renormalize the molecular orbitals. Orbitals 
with energy close to the Fermi level are pulled more to- 
wards the surface while higher and lower lying orbitals 
are pushed away. As a result, not only the energies but 
also the life-times (inverse spectral width) of the molec- 
ular resonances are affected by the xc self-energy. 

The paper is organized as follows. In Sec. |TI]we briefly 
review the general quasiparticle theory including a trans- 
parent definition of the QP states and a simple proof 
of the equivalence between this definition and the QP 
equation. In Sec. IIIII we introduce the metal-molecule 
interface model and in Sec. IIVI we explain the method 
used to calculate the QP energies and wave functions. In 
Sec. [V] we present the results and discuss implications 
for modeling of charge transport in molecular junctions. 



II. QUASIPARTICLE THEORY 

In this section we review the concept of the QP wave 
function and discuss its physical meaning. For simplic- 
ity we shall make the assumption that the system under 
consideration is finite and the relevant excitations are 
discrete^. 

We denote the A^-particle groundstate and the excited 
states by \^q) and \^f^), respectively. The occupied and 
unoccupied QP orbitals are denoted ) and \ipt), re- 
spectively. The QP orbitals belong to the single-particle 
Hilbert space and are defined through their matrix ele- 



ments with a general orbital \<f>): 

{<t>\%Y = (2) 

Wt) = (3) 

where and ct annihilates and creates an electron in 
the orbital \<j>), respectively. The real space representa- 
tion of the QP wave functions are obtained by setting 
\<ft) = \r) in the above equations. The QP wave functions 
defined above are also sometimes referred to as Lehman 
amplitudes or Dyson orbitals. 
The QP energies are defined by 

e~ = E» (4) 
4 = E^ 1 E» (5) 

They represent the excitation energies of the N ± 1 
relative to Eq and thus correspond to electron addi- 
tion/removal energies. 

The definition of the QP wave functions given in Eqs. 
(|2I3[) is not very transparent at first sight. A more trans- 
parent definition of the QP states can be obtained by 
noting that the projection 

\(*Z + %\*$)\ 2 ... 

is maximized exactly when \<p) equals \tp^). In other 
words, \ip£) is the orbital that makes c^I^q) the best 
approximation to the excited state Similarly, 
is the orbital that makes c^I^q) the best approx- 
imation to the excited state l^f^" 1 ). Consequently, the 
QP wave function is the single-particle orbital that best 
describes the state of the "extra" electron/hole in the 
excited state |4'^ r±1 ). In general, the QP states are non- 
orthogonal and their norm lies between and 1. The 
norm is a measure of how well the excited many-body 
state can be described as a single-paticle excitation from 
the groundstate. 

In the special case of non-interacting electrons, the QP 
wave functions have norms exactly 1 or 0. The former 
correspond to excitations where one extra particle has 
been added to the groundstate Slater determinant. In 
this case the QP wave functions coincide with the nor- 
malized eigenstates of the one-electron Hamiltonian. The 
QP states with zero norm correspond to all other types 
of excitations. 

It should be noted that the term "quasiparticle state" 
is often used only for those \tpp ) whose norm is close 
to 1 while other states (those corresponding to collective 
excitations) are referred to as " satellites" . In the present 
work we shall only consider QP states with norms very 
close to 1. 

In the case where E^ 1 belongs to the discrete spec- 
trum of the many-body Hamiltonian, it can be shown 
that tp^; and are solutions to the QP equation (JTJ) . 
In this case the norm of the QP state is given by Z = 



3 



(1 — dT, xc (e^)/de) l . The definition of QP states belong- 
ing to the continuum is a bit more tricky 28 , but this has 
no consequences for the present work. 



III. MODEL 
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FIG. 2: (Color online) a) Lattice model of a molecular wire 
interacting with a metal surface via hopping and Coulomb 
interaction, b) The QP density of states (spectral function) 
of the adsorbed molecule. The different broadening of the 
levels reflects the difference in the shape of the QP orbitals. 

We consider a four-unit paraphenylene molecular wire 
connected to a metal surface, see Fig. [5] The parapheny- 
lene molecule is described by a Pariser-Parr-Pople (PPP) 
Hamiltonian 2 ^ 



E 



J(J' ^1<J 1 



(7) 



where q CT (ci CT ) creates (annihilates) and electron at site i 
(p z orbital of carbon atom i) with spin a . The first term 
describes nearest neighbor hopping of strength t = —2.4 
cV. In the second term, Vy is the long range Coulomb 
interaction acting between all sites on the molecule and 
for which we use Ohno's parametrization 26 



V iS = 
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(14.4/C/) 2 + K? 



(8) 



where Rij is the distance between atoms i and j (in A) 
and U = 11.26 (in eV). We note that the PPP model with 
the parameters used here in general provides an accurate 
description of the low-lying excitations in ir conjugated 
systems 2 ^. 

As a qualitative model of the metal surface we use a 
semi-infinite one dimensional tight binding lattice. We 
use a large hopping parameter of to = —5.0 eV between 
the sites of the chain to simulate a broad featureless band 
and we set Ep = corresponding to half-filling. The 
last site on the metallic chain is coupled to the nearest 
carbon atom of the molecule by the hopping parameter 
thyb = 1-0 eV (see inset of Fig. @|. Coulomb interactions, 
Vij, as defined in Eq. ([5} are included between the last 
site on the chain and all the sites of the molecule. We 



set the distance between the last site of the chain and 
the contacting carbon atom to 3 A. With these parame- 
ters, the model yields realistic image charge shifts in the 
range 0.2-1.0 eV depending on the spatial form of the 
orbita l 18 ' 19 . We stress that the Fermi energy lies in the 
middle of the gap between the highest occupied molecular 
orbital (HOMO) and lowest unoccupied molecular orbital 
(LUMO) such that the contacted molecule remains in a 
closed-shell configuration far from the strong correlation 
Kondo regime. 

We note that the use of a ID chain to simulate the 
metal surface is clearly not adequate for quantitative 
computations. In particular, it cannot be used to de- 
scribe the reorganization of electrons in the metal surface. 
However, from the viewpoint of the molecule it captures, 
at least qualitatively, all the aspects of the image charge 
effect in a real metal-molecule junction. In particular, 
the effect that a charge added to the molecule induces an 
image charge in the metal (change in the occupation of 
the last site of the ID chain) which acts back on the elec- 
trons of the molecule. We also note that Delta-SCF as 
well as ab-initio GW calculations for a benzene-diamine 
molecule adsorbed on an adatom, or small pyramid tip 
structure, on a gold surface, have shown that the image 
charge induced in the metal is largely confined to the 
adatom and is thus quite localized, see Figure 1 in Ref. 

M 



IV. METHOD 

To obtain the QP wave functions and energies we cal- 
culate the single-particle Green function following the 
method described in Ref. H(J Briefly, the Green func- 
tion of the contacted molecule is calculated from 



Gij(0j) = [U-H - EhybM - SxcM] 



(9) 



where Hq is the non-interacting part of the molecular 
Hamiltonian including the Hartree field and Shyb is an 
embedding self-energy accounting for the coupling to 
the semi-infinite chain. In this work the xc self-energy 
is evaluated using either the Hartree-Fock or the GW 
approximation. Unless explicitly stated, the GW self- 
energy is evaluated fully self-consistently. The energy 
dependence of G and Sgw is sampled on a uniform grid 
u n = En + iil where 77 = 0.01 is an imaginary infinitesimal 
and e n ranges from (in eV) -100 to 100 with a spacing of 
rj/2 = 0.005. 

We have previously shown that the GW approxima- 
tion yields QP energies of molecules described by PPP 
models in good agreement with exact diagonalization re- 
sults with an average deviation of the lowest QP energies 
of less than 5%2£. In that work we also showed, us- 
ing a measure for the degree of correlation based on the 
entropy of the reduced density matrix, that PPP mod- 
els are significantly less correlated than Hubbard models 
with the same interaction strengths (obtained by remov- 
ing all long-range interactions Vij with i ^ j from the 
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PPP model), explaining earlier studies which concluded 
that GW does not perform well for Hubbard cluster s 31 ' 32 . 

The Green function is related to the QP states and en- 
ergies via its Lehmann representation 5 . Using this rep- 
resentation, the spectral function, A(lu) = (i/ir)[G(u) + 
irj) — G{lo — irj)], projected onto the sites of the 

molecule can be written as 



A zj (uj) 



se{+,-} 



(10) 



We identify the molecular QP energies, e" 10 ', of the 
molecule as the peaks in A(uj). The corresponding QP 
orbital (precisely, the projection of the QP orbital of the 
infinite metal-molecule system onto the molecule) is ob- 
tained as the unique solution to the eigenvalue equation 



E 

j g mol 



A ij -(C° 1 )01C oi ) = A{i|< ol > 



mol \ 



(ii) 



with A / 0. In other words, the QP orbitals are ob- 
tained by diagonalizing the matrix A at its peak ener- 
gies, and picking the eigenvector corresponding to the 
largest eigenvalue (in practice we find the largest eigen- 
value is about 10 3 times larger than the second largest). 
We stress that this method of obtaining the QP wave 
functions is equivalent to solving the QP equation ([1}, 
but is more convenient for systems with open boundaries. 
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FIG. 3: (Color online) (a) The polyparaphenylene molecu- 
lar wire, (b) Band structure of the infinite polyparapheny- 
lene molecular wire calculated with DFT(LDA) and the PPP 
model Hamiltonian with interactions treated at the Hartree 
level, respectively. The two highest valence bands (A and B) 
and two lowest conduction bands (A* and B*) are indicated. 



V. RESULTS 

In Fig. [3] we show the band structure of the infinite 
poly-paraphenylene wire obtained from a DFT-LDA cal- 
culation (red) and the PPP model with interactions de- 
scribed at the Hartree level (blue) . We have verified that 
the LDA xc-potential merely provides a constant shift 
of all the bands, and thus the two levels of approxima- 
tion are directly comparable. We conclude that the PPP 
model yields a reliable description of the it bands of poly- 
paraphenylene. The bands denoted A and A* are mainly 
composed of the HOMO and LUMO orbitals of the ben- 
zene units, while the narrow B and B* bands are formed 
by the HOMO-1 and LUMO+1 benzene orbitals, respec- 
tively. 

In Fig. 2]we show the projected density of states of the 
4-unit paraphenylene molecule coupled to the metallic 
chain 



PDOS(e) - J2 A «( £ )- 



(12) 



iGmol 



For five molecular levels we indicate the shift due to cor- 
relations (mainly the image charge effect) by horizontal 
arrows. The single-particle orbitals obtained from HF 
are shown to the left while the QP orbitals from GW are 
shown to the right. The weight of all the depicted QP 
orbitals is very close to 1 and they are essentially orthog- 
onal indicating that the single-particle picture applies. 

The HF orbitals are completely delocalized over the 
molecule and are essentially identical to the orbitals of 
the freemolecule. This is in sharp contrast to the orbitals 
derived from GW which are localized on different parts 
of the molecule. The localization of the QP wave func- 
tions occur because of the interaction between the hole 
on the molecule and the image charge that it induces in 
the metal surface. This is a highly non-local correlation 
effect and is completely missed by the HF approximation. 

It is clear that the orbitals belonging to the narrower B 
band become more localized than the orbitals belonging 
to the wider A band. This is because it is energetically 
cheaper to redistribute the orbitals of a narrow band. Fo- 
cusing on the B1-B4 states we observe a clear trend in 
the localization: The closer the energy of an orbital is 
to Ep, the closer to the surface is the orbital localized. 
Due to the image charge effect it is always energetically 
favorable for the hole to reside closer to the surface. On 
the other hand, the QP orbitals should remain (almost) 
orthogonal, at least when the QP picture applies as is the 
case here, and this prevents that all orbitals contract to- 
wards the surface. We recall from Eq. ((4]) that occupied 
QP orbitals, ip~, lying closer to Ep correspond to many- 
body excitations, Vt^ -1 , with lower energy. The observed 
trend in the localization then follows from the variational 
principle applied to the may-body states 

The unoccupied orbitals are affected by the metal sur- 
face in a similar way with orbitals lying closer to Ep 
becoming localized more towards the surface and experi- 
encing a larger energy shift toward Ep. The fact that the 
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FIG. 4: (Color online) Projected density of states of a 4-unit paraphenylene molecular wire coupled to a metal surface (middle 
panel). The red and green curves show results obtained at the Hartree-Fock and GW level, respectively. For five energy levels 
the shift due to correlations (mainly the image charge effect) is indicated by arrows, and the corresponding orbitals are plotted 
in the left (Hartree-Fock) and right (GW) panels. Orbital Al is the HOMO and belongs to band A in Fig. |3]while the orbitals 
B1-B4 belong to the narrow band B. The orbitals were constructed by superposing p z orbitals with weights given by the discrete 
wave functions of the model. 



sign of the image charge shift of the energy of empty and 
occupied orbitals is different shows that the effect cannot 
be mimicked by a local —1/z potential. Such a poten- 
tial would shift all orbitals in the same direction (down- 
wards). To mimick the image charge potential would 
instead require a non-local potential of the form 



1/zPocc - l/zP f 



empty 



(13) 



where P occ and P cmp ty project onto the subspace of oc- 
cupied and unoccupied molecular orbitals, respectively. 
From this property of the image charge potential, it is 
clear that the effects presented in Fig. [4] cannot be cap- 
tured by a local potential - even the exact xc-potential 
of DFT. 

From Fig. @]we see that not only the QP peak positions 
but also the width of the resonances is affected by the 
image charge effect. This is particularly pronouned for 
the Al orbital which becomes significantly broadened due 
to the increased weight of the orbital at the carbon atom 
connected to the metallic chain. For the B orbitals, which 
have very little weight on the contacting carbon atom, 
the small increase in the GW peak width relative to HF 
comes from the (small) imaginary part of the GW self- 
energy. 

The QP energies include correlations in addition to 
the exchange effects described at the Hartree-Fock level. 
The correlation energy contains contributions from the 
Coulomb interactions internally on the molecule (internal 
screening) as well as the interactions between metal and 
molecule (image charge screening) . From calculations for 
the molecule in the gas-phase we have verified that the 
contribution from internal screening is almost the same 
(between 0.4 and 0.6 eV) for the different molecular or- 
bitals. Hence, apart from this constant, the difference 
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FIG. 5: (Color online) The size of the image charge effect is 
plotted for the orbitals A1-A4 and B1-B4 against the center of 
the QP orbital along the axis of the molecule, (x) . The degree 
of localization quantified as the second moment {(x — {x}) 2 ), 
is indicated by a horizontal line for each orbital. 



between the QP energy and the HF energy represents 
the shift in the energy level due to the image charge ef- 
fect. 

In Fig. [5] we plot the image charge shift for the A 
and B orbitals against the center of the QP orbital along 
the axis of the molecule (the x-axis). The center is de- 
fined as the first moment, (-0™ ol |i|i/ , ™° 1 )- The vertical 
dashed line indicates the center of the molecule which 
coincide with the center of the HF orbitals. For each or- 
bital, the degree of localization, quantified as the second 
moment (ij)™ \(x — (x)) 2 \ip™ ), is indicated by a horizon- 
tal line. The vertical dashed line indictaes the center of 
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FIG. 6: (Color online) The image charge shift for the A and 
B orbitals calculated using three different methods for calcu- 
lating the GW self-energy and solving the QP equation (see 
text). 



VI. CONCLUSION 

Wc have discussed the general mathematical and phys- 
ical meaning of quasiparticle (QP) wave functions in in- 
homogeneous systems. For the specific case of a molecule 
adsorbed on a metal surface, we found that the QP states 
can differ qualitatively from the orbitals obtained within 
standard independent-particle approximation as exam- 
plified by the Hartree-Fock approximation. Using the 
GW method, it was shown that image charge interac- 
tions pull the QP frontier molecular orbitals towards the 
surface. In contrast, the Hartree-Fock single-particle or- 
bitals remain delocalized and identical to those of the iso- 
lated molecule. These results are of importance for the 
modelling of energy level alignment and electron trans- 
port across metal-molecule interfaces, and should be ob- 
servable by low-temperature scanning probe experiments 
on molecules on insulating substrates. 



VII. ACKNOWLEDGEMENT 



the molecule which coincide with the center As expected 
there is a clear correlation between the size of the image 
charge shift and the orbital center, in particular for the 
highly localized B orbitals. The orbitals A2-A4 are all 
pushed slightly away from the surface. 

Fig. [5] compares the image charge shifts obtained us- 
ing three different strategies for calculating the GW self- 
energy and solving the QP equation: (i) full solution of 
the QP equation with a self-consistent GW self-energy 
(ii) full solution of the QP equation with a one-shot 
Go Wo self-energy with Go from Hartree-Fock, and (iii) 
first-order perturbation theory applied to the GoWo self- 
energy, i.e. e% p = + (^ F |E GoW , (s" F )|^ F )- As 
expected, the first-order approximation (iii) does not per- 
form well in cases where the QP orbitals deviate signif- 
icantly from the HF orbitals, i.e. for B1-B4 and Al. 
In particular, for the B1-B4 orbitals the first-order ap- 
proach predicts similar image charge shifts whereas the 
shifts obtained with methods (i) and (ii) vary due to the 
variation in the distance of the QP orbitals to the sur- 
face. As a general trend, the self-consistent treatment 
of the GW self-energy leads to larger image charge shifts 
(smaller HOMO-LUMO gaps) in agreement with findings 
for isolated molecules 3 ^. 

At first sight it might seem surprising that the image 
charge effect, which is essentially electrostatic in nature, 
is not captured by meanfield theories such as HF. The 
reason is that the effective potential defining the mean 
field Hamiltonian does not "know" how an additional 
particle on the molecule becomes screened by the metal. 
Under some limiting conditions, however, the effect of 
image charge interaction can be simulated using mean 
field methods to compute the total energy with an extra 
particle explicitly present on the molecul e 30 ! 34 . 
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